Universality Classes in Constrained Crack Growth 
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Based on an extension of the fiber bundle model we investigate numerically the motion of the 
crack front through a weak plane separating a soft and an infinitely stiff block. We find that there 
are two regimes. At large scales the motion is consistent with the pinned elastic line model and we 
find a roughness exponent equal to 0.39 ± 0.04 characterizing it. At smaller scales, coalescence of 
holes dominates the motion rather than pinning of the crack front. We find a roughness exponent 
consistent with 2/3, which is the gradient percolation value. The length of the crack front is fractal 
on these smaller scales. We determine its fractal dimension to be 1.77 ± 0.02, consistent with the 
hull of percolation clusters, 7/4. This suggests that the crack front is described by two universality 
classes: on large scales, the pinned elastic line one and on small scales, the percolation universality 
class. 
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Schmittbuhl and Mai0y [l[ were the first to study ex- 
perimentally the roughness of a crack front moving along 
a weak plane in a material loaded under mode I condi- 
tions. By sintering two sandblasted plexiglass plates to- 
gether and then plying them apart from one edge, they 
were able to follow the motion of the crack front mov- 
ing along the sintered boundary between the two plates. 
The rough crack front turned out to be self affine, i.e., its 
height-height correlation function ((h(x + Ax) — h(x)) 2 ) 
scaled as {Ax]'' where C is the rougness exponent. h(x) 
is the position of the crack front with respect to a base 
line orthogonal to the average crack growth direction and 
x is the coordinate along this base line. The roughness 
exponent was found to be C = 0.55 ± 0.05. A couple of 
years before this work, Schmittbuhl et al. [2j studied nu- 
merically a model of such constrained crack growth based 
on regarding the motion of the crack front as that of a 
pinned elastic line, the fluctuating line model. This work 
was based on an earlier idea by Bouchaud et al. Q . The 
conclusion of Schmittbuhl et al. was that the front should 
be rough and self affine, characterized by a roughness ex- 
ponent £ = 0.35 ± 0.05. This value was later refined 
to 0.388 ± 0.002 by Rosso and Krauth [j. The large dis- 
crepancy between the numerical and experimental results 
- the latter having been refined to ( = 0.63 ± 0.03 by 
improving the statistical analysis [H — spurred a lively 
quest for an explanation that only today seems to con- 
verge towards a satisfactory understanding of the under- 
lying physics , see e.g. [y, [Z| for a review. 

An attempt at finding an explanation for the large 
roughness exponent seen in the experiments that was not 
based on some variation on the fluctuating line model was 
forwarded by Schmittbuhl et al. [H . The underlying idea 
here was that the crack front did not advance due to a 
competition between effective elastic forces and pinning 
forces at the front, but by coalescence of damage in front 
of the crack with the advancing crack itself. Such an 
idea had been put forwards in a more general context 



by Bouchaud et al. [9( the year before. By using a grid 
of linearly elastic fibers connected to a soft elastic half 
space, Schmittbuhl et al. found a rougness exponent of 
C = 0.60±0.05, which was consistent with the experimen- 
tal results. The crack advance mechanism in this model 
was analyzed in terms of a correlated gradient percolation 
process where coalescence is the dominating mechanism. 

Recently, Santucci et al. [l(| reanalyzed data from 
a number of earlier studies, including [5J, finding that 
the crack front has two scaling regimes: one small- 
scale regime described by a roughness exponent = 
0.60±0.05 and a large-scale regime described by a rough- 
ness exponent = 0.35 ± 0.05. 

We suggest in this letter that there are indeed two 
competing mechanisms involved in generating the scal- 
ing properties seen in the roughness of the crack front: 
on small scales coalescence dominates, whereas on large 
scales, the fluctuating line picture is correct. There 
is a crossover between the regimes where either of the 
two mechanisms dominate associated with a well-defined 
crossover length scale. 

Our numerical work strongly suggest that the coales- 
cence mechanism seen at small scales is controlled by the 
ordinary percolation critical point, leading to percolation 
exponents. This is in constrast to Schmittbuhl et al. Q 
who suggested a new percolation-like universality class, 
but consistent with recent theoretical work on fracture in 



the fuse model [11 1 . 

We note an earlier attempt at explaining the two scal- 
ing regimes seen in Ref. [10( based solely on the fluctu- 
ating line model 12 1. Here Laurson et al. 12| relate the 
crossover to the Larkin length scale of the crack front 
[l3| . This is in this context the length scale at which the 
roughness of the front is comparable to the correlation 
length inherent to the pinning disorder. 

Our numerical system is a fiber bundle model [13] 
closely related to the one used in [§] and introduced by 
Batrouni et al. the year before 15|. L x L elastic fibers 
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are placed in a square lattice between two clamps. One 
of the clamps is infinely stiff whereas the other has a fi- 
nite Young modulus E and a Poisson ratio v. All fibers 
are equally long and have the same elastic constant k. 
We measure the position of the stiff clamp with respect 
its position when all fibers carry zero force, D. The force 
carried by the fiber at position where i and j are co- 
ordinates in a cartesian coordinate system oriented along 
the edges of the system, is then 



f{i,3) = - fc ( u (*j) _ D ) 



(1) 



where Uuj\ is the fiber's elongation. The fibers redis- 
tribute the forces they carry through the response of the 
clamp with finite elastisity. The redistribution of forces 
is accomplished by using the Green function connecting 
the force ff m ,n) acting on the clamp from fiber (m, n) 
with the deformation at fiber [IH 
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where a is the distance between neighboring fibers. 
ruj) — T7 TO ,„) is the distance between fibers and 
(m, n) and the integration runs over the a x a square 
around fiber (m,n). This equation set is solved using a 



Fourier accelerated conjugate gradient method [13 

We note that the Green function, Eq. J2bj, is propor- 
tional to (Ea) -1 . The elastic constant of the fibers, k, 
must be proportional to a 2 . The linear size of the sys- 
tem is aL. Hence, by changing the linear size of the 
system without changing the discretization a, we change 
L — > XL but leave (Ea) and k unchanged. If we on 
the other hand change the discretization without chang- 
ing the linear size of the system, we simultaneously set 
L -> AL, (Ea) -> X(Ea) and k -> k/X 2 . We define the 
scaled Young modulus e = (Ea)/L. Hence, changing e 
without changing k is equivalent to changing L — and 
hence the linear size of the system — while keeping the 
elastic properties of the system constant [l9| ■ 

Bonds are broken by using the quasistatic approach 
[20l | . That is, we assign to each fiber a threshold 

value t(i j) . Bonds are broken one by one by each time 
identifying m,axuj-\ (fn j) jti $ j\ ) for D = 1. This ratio is 
then used to read off the value D at which the next fiber 
breaks. 

In the constrained crack growth experiments of 
Schmittbuhl and Mal0y [l(, the two sintered plexiglass 
plates were plied apart from one edge. In the numer- 
ical modeling of Schmittbuhl et al. 8], an asymmetric 
loading was accomplished by introducing a linear gra- 
dient in D. Rather than implementing an asymmetric 



loading, we introduce a gradient in the threshold distri- 
bution, tuj) = gj + i"(i t j) , where g is the gradient and 
Taj) is a random number drawn from a flat distribution 
on the unit interval. In the limit of large Young mod- 
ulus E. this syshnn becomes equivalent to the gradient 



percolation problem [21 



In order to follow the crack front as the breakdown 
process develops, we implement the "conveyor belt" tech- 
nique 23, 24| where a new upper row of intact fibers is 



added and a lower row of broken fibers removed from the 
system at regular intervals. This makes it possible to 
follow the advancing crack front indefinitely. The imple- 
mentation has been described in detail in 
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Figure [T] shows two examples of typical crack fronts 
representative of a stiff (high e = 0.8) and a soft system 
(low e = 2xlCP 3 ). In both cases, the front is viewed 
from above (and propagation corresponds to an upwards 
movement in the figure). The two fronts are quite dif- 
ferent. Even though we change the Young modulus E 
between the two crack fronts, this is equivalent to chang- 
ing the scale at which the fronts are viewed: Small e 
corresponds to a large scale and vice versa. 

We identify the crack front by first eliminating all is- 
lands of surviving fibers behind it and all islands of failed 
fibers in front of it. We measure "time" n in terms of 
the number of failed fibers. After an initial period, the 
system settles into a steady state. We then record the 
position of the crack front j = j(i,n = 0) after having 
set n — 0. We then define the position at later times 
n > relative to this initial position, 



h i( n ) = j{h n ) -j(i,n= 0) . 



(3) 



This is the same definition as was used by Schmittbuhl 
and Mal0y [l|. The front as it has now been defined will 
contain overhangs. That is, there may be multiple values 
of hi(n) for the same i and n values. We only keep the 
largest ft,*(n), i.e. we implement the Solid-on-Solid (SOS) 
front. 

We define the average position of the front as 
(h(n)} = Yli=i hi(n)/L and the front width as w(n) 2 = 
Eti(hi(n)-(h(n))) 2 /L. 

We will in the following explore the model going from 
large to small values of the Young modulus E while keep- 
ing k fixed. As already discussed, this is equivalent to 
changing the system size L while keeping E fixed. 

We analyze the fracture fronts in the following using 
the average wavelet coefficient (AWC) method [25L l26j|. 
We transform the front hi(n) using the Daubechies-4 
wavelet, hi(n) — > w a ,b where a is scale and b is posi- 
tion (and we have suppressed the n dependency). We 
then average \w a ,b\ for each length scale a over position 
b, and if hi(n) is self affine with a roughness exponent £, 
we have 



W(a) = <K b |) ~a< +1 / 2 



(4) 




FIG. 1: Parts of the crack fronts obtained in the simulations. The color black represents broken bonds and white represents 
unbroken bonds. To the left, a) is an example from a system with a high rescaled Young modulus e being driven forwards 
(upwards) primarily by coalescence with damage forming ahead of the crack front, (b) is an example from a system with a low 
e. This system is being driven forward by damage forming on the crack front. Small e = 2 x 10~ 3 is equivalent to a larger 
system than one with a large e = 0.8. 




FIG. 2: Averaged wavelet coefficient W(a) based on trans- 
forming hi vs. a for e = 7.8125 x 10~ 4 , L = 256 (g = 0.00625, 
averaged over 4200 fronts) and L = 512 (g = 0.003125, av- 
eraged over 850 fronts). The slope of the straight line is 
0.39 + 1/2. 



We start by considering systems with small scaled 
Young modulus e. This is a soft system at small length 
scales — or equivalently, a stiffer system at large length 
scales. We set e — 7.8125 x 10~ 4 . The fronts in this 
regime have an appearance as in Fig. [Tp. In Fig. [2j we 
plot the averaged wavelet coefficient W(a) against the 
scale a. The data follow a power law, W(a) ~ a °- 39 + 1 / 2 ; 
leading to a roughness exponent of 



C4 



0.39 ±0.04 



(5) 



entirely consistent with the large scale rougness exponent 



measured by Santucci et al. [lfj, C+ = 0.35 ± 0.05. 

We now turn to large scaled Young modulus, i.e. stiff 
systems — or equivalently, softer systems on small length 
scales. Hence, fronts appear as in Fig. [T^,. We set in the 
following e = 3.125. The corresponding plot of averaged 
wavelet coefficient W(a) vs. scale a is shown in the insert 
in Fig. [3^- The different curves correspond to different 
gradients a and L. Repeating the analysis of Sapoval et 
al. 2lj, l27[ for gradient percolation, we assume that the 



it. In the case of gradient percolation, this correlation 
length is related to the gradient through the relation £ ~ 
g-v/i^+v) = g -4 / 7 , where v = 4/3 is the percolation 
correlation length exponent. As was argued in Ref. 
the same analysis may be repeated for correlated systems 
such as the present one, but with a possibly different 
v. We would then expect data collapse in Fig. [3ji by 
rescaling W(a) -> W(a)/£ = W(a)/g~ a and a ->■ a/£ = 
a/g^ 13 , and where a = (3 = u/(l + v). We show in Fig. 
[3ja,, the data collapse that ensues when choosing a — 
P = 4/7, i.e., the ordinary uncorrelated percolation value 
v/(l + is), where v = 4/3. There are two distinct regions 
in the figure. For large values of a/g~^, W(a)/g~ a is 
independent of aj ' g~^ '. Hence, the front has the character 
of uncorrelated noise. We will discuss this further on. For 
small a/g~P , the data follow a power law. We show in the 
figure (Fig. [5^) a straight line with slope 7/6 = 2/3 + 1/2. 
Hence, the data are consistent with the fronts being self 
affinc on these scales with roughness exponent 



C- = 2/3 



(6) 



a value consistent with gradient percolation 27[ — and 
consistent with the experimental value reported by San- 
tucci et al. 0, C- = 0.60 ± 0.05. 

In order to further the analysis, we follow the proce- 
dure in 2jj by smoothening the fronts by removing the 



jumps due to the overhangs through the transformation 

hi(n) -+ h^ln) = 

i 

^ sgn[/i m+ i(n) - h m (n)] \h m+1 (n) - h m (n)\ k , 



m— 



(7) 



front has an isotropic correlation length £ associated with 



when k — > [28(. Hence, all steps are equal to one in 
h®(n), all overhangs have been removed. We ensure by 
this procedure that the scaling properties of the rough- 
ness are due to self affinity and not due to the overhangs 
which are prevalent in this regime. 

We plot in Fig. [3p the rescaled W(a)/g~ a based on 
transforming h® vs. the rescaled a/g~@. We have set 
a = (3 = 4/7, consistent with ordinary gradient percola- 
tion. The two straight lines that have been added to the 
figure have slopes 2/3 + 1/2 = 7/6 and 1/2 + 1/2 = 1 
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FIG. 3: a) Averaged wavelet coefficients W(a) of the SOS 
front hi vs. scale a. The straight line has slope 2/3 + 1/2 = 
6/7. The insert shows the unsealed data. The data are based 
on e = 3.125, L = 256 (g = 0.004, 0.005, 0.01, 0.02 and 0.03). 
b) Average wavelet coefficients W(a) of the SOS fronts h® 
where the overhangs have been removed. The long-dashed 
line has slope 2/3 + 1/2 = 6/7 and the short-dashed line has 
slope 1/2+1/2 = 1. In both a) andb), a = p = 4/7. The data 
are based on e = 3.125, L = 64 (g — 0.018 and 0.02, averaged 
over 1426 fronts), L = 128 (g = 0.008, 0.01, 0.03, and 0.05, 
averaged over 100 to 500 fronts), L = 256 (g = 0.004, 0.005, 
0.008, 0.01, 0.015, 0.02, 0.025, 0.03, and 0.035, averaged over 
30 to 500 fronts), and L = 512 (g = 0.002, averaged over 50 
fronts) . 



respectively. On small scales, the fronts are then self 
afnne with a roughness exponent = 2/3 — the gra- 
dient percolation value [27j |. On larger scales and with 
the overhangs removed, one would naively have expected 
to observe the fluctuating line regime characterized by 
a roughness exponent C+ = 0.39. However, by remov- 
ing the overhangs, the effective roughness exponent one 
measures is max(l/2,£+), which in this case is 1/2 (28j |. 
Fig. [3Jd, then, shows the crossover from a roughness ex- 
ponent consistent with ordinary gradient percolation to 
a plain random walk exponent which is a result of the 



smoothening process. 

Hence, we conclude that the scaling properties seen 
on small scales are consistent with uncorrelated gradient 
percolation. The analysis of Schmittbuhl et al. [1| sug- 
gested a correlated gradient percolation process. This 
is one step further. We suggest that the process is an 
uncorrelated gradient percolation process on small scales 
where coalescence is the dominating process. 

Roughness exponents are notoriously difficult to mea- 
sure. The data presented in Fig. [3] are not of sufficient 
quality to warrant firm conclusions on the small-scale 
universality class by themselves. We therefore in the 
following measure the fractal dimension of the fronts. 
Leaving the SOS assumption, we now follow the front as 
shown in Fig. [TJ The front has a length I. For small val- 
ues of the scaled Young modulus e, there are no (or very 
few) overhangs and we expect I to be proportional to L: 
it is not fractal. However, for large e where overhangs 
are prevalent, we do expect it to be fractal. We assume 
that there is a correlation length £ and that the front is 
fractal up to this scale. Hence, the length of the front I 
then scales as 



I 



L 



(8) 



From Sapoval et al. [2l[, w e know that £ ~ g -"/Q+ v ). 
We now set g = c/L, where c is a constant. Hence, we 
find 



(9) 



If we are dealing with ordinary gradient percolation, v = 
4/3 and D f = 7/4 Hence, we expect (vD f + l)/(i/+ 
1) = 10/7 « 1.43. We show in Fig. gj l/L w / 7 as a 
function of the scaled Young modulus, e. For large values 
of e, there is excellent data collapse. For small values of 
e, there is data collapse when l/L is plotted against e, 
see the insert of Fig. HI indicating that the front is not 
fractal in this regime. We are seeing the crossover from 
the coalescence-dominated regime to the fluctuating line 
regime as e moves from large to smaller values. Since 
e = E/L, we may keep E fixed and change e by changing 
L. Hence, the coalescence regime is a small-scale regime 
whereas the fluctating line regime is a large-scale regime. 

If we work backwards and search for the best scal- 
ing exponent to produce data collapse in Fig. [4] we find 
l/L 1AA . This gives Df = 1.77 ± 0.02, entirely consistent 
with the percolation value 7/4. 

Hence, we have through the use of a single model based 
on the fiber bundle model with a soft clamp and a gradi- 
ent in the breaking threshold, been able to identify two 
mechanisms by which the fracture front propagates. On 
small scales, it is coalescence of damage the dominates, 
whereas on large scales the front advances as descibed by 
the fluctuating line model. The coalescence regime is in 
the universality class of ordinary percolation. 
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FIG. 4: Length of the fracture front I scaled by system size 
plotted against effective Young modulus e. The main figure 
shows data collapse for l/L 10 ^ 7 for large e values and the 
insert shows data collapse for small e values for l/L . We set 
g — 1.6/ L. The data are based on 425000 fronts for L = 32, 
286482 fronts for L = 64, 282022 fronts for L = 128 and 32117 
fronts for L = 256. 
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